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(54) Method of stabiUng spectra In spectometry 

(57) A method of quarMatove analysts of a specific 
component of an abf«ct to t>e measured comprises 
steps of measuring an mnmgf spectrum of light trans- 
mitted through or reflected t»om tt>e object, dividing the 
energy spectrum into a piuraifty of wavelength domains, 
thereby obtaining a pkf alty of partial energy spectra, 
normalizing the plurality of partal energy spectra within 
each wavelength doman usmg an energy measured at 
a predetermined wavelength contained in each wave- 
length domain, and pertgrng^g me quantative analysis 
by multivariate analyse usmg me plurality of partial 
energy spectra having been normalized. 
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Description 

BACKGROUND OF THE INVENTION 

5 1 . Reltf of the Invention 

The present invention relates to methods of stabilizing spectra in spectrometry and more particularly methods of 
stabilizing measurements by reducing in numerical computation the fluctuations of the spectral baseline due to changes 
in temperature, source voltage, and the like within the spectroscope and the fluctuations of transmittance due to 
10 changes in the incidence angle of the light applied to the measured object. 

2. Description of the Related Art 

In general, quantitative analysis for the concentration of a specified component can use an analytical curve 
IS obtained beforehand by measuring the energy spectrum of light transmitted through or alterratively reflected by a sam- 
ple whose concentration is already known. As methods of obtaining the analytical curve, generally known are a method 
in which a single wavelength is used and a method in which two or nnore wavelengths are used. 

However, fluctuations of the spectral baseline occur in the spectrum caused by changes in temperature inside and 
outside the spectroscope, changes in the voltage of the power supply, and the like. These fluctuations of the baseline 
20 reduce the accuracy of the quantitative analysis, so that it is necessary to eliminate the error components in the spec- 
trum. 

As methods of reducing the baseline fluctuations by numerical computation, dout)l&wavelength baseline compen* 
sation, digital differentiation, and the low-frequency filtering Fourier transform are generally known. Each of these meth- 
ods is briefly described in the following. 

25 

DOUBLE-WAVELENGTH BASELINE COMPENSATION. 

The method of double*wavelength baseline conrpensation uses two wavelengths and X2 at which the absort>- 
ances are constant for any solution composition to cancel the t>aseline component out. Specifically, the two wave- 
30 lengths and X2 respectively satisfy the following equations (1 ) and (2). 



^a/ x/C/'^constl (Constant), (1) 

35 

where 

ai,xi • at)sorptivity at wavelength 
q: concentration of component i, 
40 n: the number of components in the sdutton, 

/: cell length. 



Z**/A2^/ = const2 (Constant). (2) 
45 /-I * * 



where 

• aj.x2' ^^^n^tivity at wavelength Xgt 
50 q: corx:entration of component /, 

n: the number of components in the solution, 
/: celt length. 

Therefore, if the absort)ances of a measured absorption spectrum f(X) at wavelengths and X2 are respectively 
constS and oonst4, then the compensated spectrum is obtained by subtracting 
55 (const4 - const3)/(X2- X^)x(x -X^)-^ constS from fix) . 
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DIGfTAL DIFFERENTTIATION. 

The digital differentiation method reduces or eliminates a drift at a low frequency by digitally differentiating absorb- 
ance with respect to a wavelength. The order of differentiation is generally first or second. The constant component 
5 independent of wavelength can be eliminated by the first order differentiation. The first-order drift with respect to the 
wavelength can be eliminated by the second order differentiation. In digital differentiation. absort)ance is differentiated 
with respect to the wavelength, so that differentiation higher than the second order greatly distorts the signal. Therefore, 
the digital differentiation of the order higher than third is rarely used. 

10 LOW-FREQUENCY FILTERING FOURIER TRANSFORM. 

The method of the low-frequency filtering Fourier transform that Fourier transforms the spectrum, i.a, converts the 
energy spectrum from the time domain to the frequency domain, to filter out the low-frequency components that are the 
causes off the baseline shift. 

15 The method of double-wavelength compensation requires the two wavelengths and X2 ^ compensation 
to respectively satisfy the above equations (1) and (2). However, in nuilti-component solutions of more than two com- 
ponents, there often do not exist the two wavelengths that satisfy the above equations (1) and (2). In these cases the 
method of dout)le-waveiength compensation cannot be used. 

The method of digital differentiation enlarges high-frequency noise such as random noise by digital differentiation. 
20 Therefore, the bandwidth of the baseline fluctuations has to be wider than the bandwidth of the light absorbed by sub- 
stances. Further, the derivative spectra are distorted by digital differentiation. 

In the method of the low-frequency filtering Fourier transform, there is a problem that the information about light 
absorption by substances in a low-frequency domain is lost by filtering together with kjaseline fluctuation components. 
Further, there is a problem common to digital differentiation and the low-frequency filtering Fourier transform. That 
25 is. if the energy F is a function of two variables, time t and wavelength X. then the total differential of E(X. i) is expressed 
by the following equation (3). 

dE = (dEldfK),(fK + {dEldO ^dt (3) 

30 The prior methods of digital differentiation and the low-frequency filtering Fourier transform are both operations on 
wavelength X. Therefore, the second term (dEldt)x of the equation (3), which is a function of f. can not be reduced or 
eliminated. 

SUMMARY QP TtiE INV^ffnON 

35 

The object of the present invention is therefore to provide methods of stabifizing spectra in spectrometry such that 
spectral error components deperxiing on time are reduced, and fluctuations of transmittance due to changes in the inci- 
dence angle of the light applied to the measured object are also reduced. 

The present invention is based on the following consideration. Errors in an energy spectrum due to disturt^uices 
40 such as changes in temperature of the spectroscope, changes in the voltage of the power supply, and the like can gen- 
erally be expressed as a drift, that Is. a temporal change of the spectrum. Therefbre, the intensity of transmitted light 
/(X. t) at wavelength X and time t and the intensity of the absorbed light at wavelength Xr and time t can be respectively 
expressed by the following equations (4) and (5). 

45 

/(X.O = k{XJ) . /o(X.f 0) • p{X.e) • exp(.5^a,(X)c,-/l, (4) 



so l{X,J) - k{X,.t) . lo(X,.to) 'P{^r^)'^l'^^ii>^r)cA' (5) 

# 



Here. 

55 /(X. t): intensity of transmitted light at wavelength X and time t, 

/(Xp f): intensity of transniitted light at wavelength X, and time /. 
ib(^. h)' intensity of incident light at wavelength X and time to, 
h(K h)' ir^ensity of incident light at wavelength X, and time ^0. 
p(X. e): surface transmittance of light incident at angle 6 at wavelength X. 
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pikf, 0): surface transmittance of light incident at angle 0 at wavelength kp 
ajik). absorption coefficient of corrponent / at wavelength X. 
a,(Xf): ak)sorption coefficient of component / at wavelength Xp 
6: incidence angle of light with the surface of the measured object 
q: concentration of component /, and 
/: effectrve path length. 

Further. A(X. 0 the equation (4) is the ratio of the intensity of incident light at wavelength X and time f to the inten- 
sity loiX, to) of inadent light at wavelength X and time ^ Therefore. 

k{X,to)^^^ (6) 

Similarly. k(K^ f) nihe equation (5) is the ratio of the intensity of incident light at wavelength Xr and time / to the intensity 
lo(Xr, to) of incident bght at wavelength x^ and time Jb- Therefbre. 

The ratio of the equebon (4) to the equation (5) yields the ratio of energy intensities expressed by the following 
equation (8). 



Here, we assume 

aj(X)^a|(X,). (9) 
The total differential of * (iu r) can be expressed by the following equation (10). 

dk{x, t) = {dk/dX) ^X + (dk/dt) ^dt. (10) 

Here, (dklcfK)x expresses tht dependency on wavelength of the ratio k{X, fjxA intensities of incident ligfit, and 
{dkldt)x expresses the rate of ctiange with respect to time of the ratio k{X, t) of intensities of the incident light. 

If the ratio of ttie intensity of modent light at time / and wavelength X, i.e.. /c(X, f)^lo(X, t), to tiie intensity of incident 
light at time f and wavelengm i • . kix, Ox^(x. t), is denoted by p(X. Xf. t), then 



Now, if the wavelength donwi a divided into subintervals, and wavelengths X and x^ are limited to the fth subinter- 
val [Aj-min. Xj-max]. then the term (<>^i?x)x. which expresses ttie dependency on wavelength of the ratio k(X, t), can be 
approximated by the foUowng ( 12) 

(O/c/OX), «0. (12) 

Therefore. 

/f(X.O*MO. (13) 

Consequently. p(X. X,. t) can be expressed by tiie following equation (14). so tiiat it can t>e also expressed as the 
equation (15). 

Slf,-^ "4, 
"MO*/o(>t^o) 
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From the above equation (15), it can be seen that en-or components due to the time dependent baseline drift in a 
spectrum can be reduced or eliminated. 
£ Next, it is shown that the fluctuations of transmittance due to changes of the incidence angle of light can be reduced 
by obtaining the ratio of energy spectra in subdomains. 

Accading to Fresnel's law. the surface transmittance p{X, O) when light is incident on the measured object at an 
incidence angle e, at wavelength X and the surface transmittance p(Xp 8) when light is incident on the measured object 
at an inciderKe angle Oj at wavelength can respectively be expressed by the following equations (16). (17). (18). and 
10 (19). 

cose, / 2n.,cos0; 



15 



20 



25 



''^^•"'i^n^j^'cose j Vn2.3^cose,+ni;^cose,>' • ^ ' 

oixe) "2.x cose, . 2n,,,cose, x2 
' "i.x cose j ^^^Lxcose^+Hg ;^,cose,/ 

''i^/cose i*U.x,cose;*n,.x^cosej " ^ ' 

Hg.x, cose, ( 2n, Jo cose, \^ 
P(>^ ,.e) ±-7r^' 55ie"| * (n ,,j,^cose ,-^n2.x,cose J ' 



Here. 



p(;i.ej, : surface transmittance of light at wavelength X when horizontally polarized light is incident on the meas- 
30 ured object at the incidence angle ei. 

p(X. e)^ : surface transmittance of light at wavelength A when vertically polarized light is incident on the measured 
object at the incidence angle 6|. 

p(Xp e), : surface transmittance of light at wavelength when horizontally polarized light is incident on the meas- 
ured object at the incidence angle Oj. 
^ PiK surface transmittance of light at wavelength Xr when vertically polarized light is incident on the meas- 

ured object at the incidence angle % 

ny ^ : index of refraction of medium 1 at wavelength A.. 
n^x- i^ex of refraction of medium 2 at wavelength X. 
n^ yjf : index of refraction of medium 1 at wavelength Xp 
40 n^ yj : index of refraction of medium 2 at wavelength X„ 

e,: incidence angle with the surface of measured object, and 
6): refraction angf e with the surface of measured ot)ject. 
Further, according to Snell's law, the following equation (20) holds. 

45 sinOj s sine,. (20) 

How, the ratio of p(X. e) to p{Xp 0) is denoted by y(X. k^, 0), i.e.. 

y(x.Xpe)«p(x.e)/p(Xpe). (2i) 

so 

tf AX Xr • X issman. then y(X. X,. e) can be approximated by 

y(X.X,. e) = y(X.X,). (22) 
55 From the equations (8). (15). and (22). we obtain the following equation (23). 

= P(^.^r) -yU.^tr) • expI-Z(a/(M-a;{X,))c,/l, (23) 
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It is seen from the above equation (23) that errors of baseline fluctuations of the spectrum due to time dependence 
and fluctuations of transmittance due to changes of the incidence angle of the incident light can be reduced at the same 
time. 

The present invention is based on the consideration desabed above. The method of the present invention is a 

5 method of spectrometry that performs quantitative analysis of a specific component of a measured object by irradiating 
the measured object with light and measuring the energy spectrum of the transmitted light or alternatively the energy 
spectrum of the reflected light. In particular, the method divides the wavelengtii domain of the measured energy spec- 
trum of the transmitted light or the reflected light into a plurality of intervals. Then the method obtains, within each inter- 
val, the ratio of each measured energy value to the measured energy value at a predetermined wavelength to stabilize 

10 the measured values of tiie spectrum. In this way, the quantitative analysis of the specific component can be achieved 
with stabilized spectral values. 

Further, according to one embodiment of the present invention, tiiere is provided a method of stabilizing spectra in 
spectrometry that performs quantitative analysis of a specific component of a measured object by in-adiating the meas- 
ured object with light and measuring the energy spectrum of the transmitted light or alternatively the energy spectrum 

15 of the reflected light. In particular, the method divides the wavelength domain of the measured energy spectrum of the 
transmitted light or the reflected light into a plursdity of intervals. Then ttie method obtains, within each interval, the ratio 
of each measured energy value to the measured energy value at a predetermined wavelength to statMlize the measured 
values of the spectrum. Then the mettiod performs multivariate analysis for ihe specto'um of the ratios. In this way. ttie 
quantitative analyst of the specific component can t^e achieved with stabilized spectral values. 

so Further, in the present invention, the specific conponent of the measured object is glucose. 

A method of ttie present invention divides the wavelength domain of the measured energy spectrum of the trans* 
mitted light or the reflected light into a plurality of intervals. Then the mettiod obtains, within each interval, ttie ratio of 
each measured energy value to the measured energy value at a predetermined wavelength to stabilize the measured 
values of the spectrum. In ttiis way. en-ors of baseline fluctuations of the spectrum due to time dependence and f luctu- 

25 ations of transmittance due to the incidence angle of the incident light can be compensated at ttie same time. 

Furtiier. a method of present invention performs multivariate analysis for the spectrum obtained by the ratios of 
energy spectra in subdomains. 

In particular, the energy spedrum of transmitted light or the energy spectrum of reflected light for glucose is com- 
pensated by taking ratios in subdomains. 

30 According to ttie present invention, errors of baseline fluctuations of ttie spectrum due to time dependence and fluc- 
tuations of transmittance due to the incidence angle of the incident light can be reduced at ttie same time. Conse- 
quentty, the accuracy and reliance of quantitative analysis in spectrometry can be improved. 

Furtiier. according to the present invention, ttie accuracy and reliance of quantitative analysis can be improved by 
multivariate analysis. 

35 Still further, according to the present invention, in ttie case of glucose being the specific component, the concen- 
tration of glucose can be detected with great accuracy by applying the mettiods of the present invention, so ttiat ttie 
measurement of ttlood sugar in human kxxfy can be obtained with great accuracy and reliance. 

BRIEF DESCRIPTION OF THE DRAWINGS 

40 

These and other objects and features of ttie present invention will become clear from the following description taken 
in conjunction with the emtxxjiments thereof with reference to the accompanying drawings, in which: 

Fig. 1 shows 1 1 energy spectra of transmitted light of commercial milk measured every 5 minutes. 
45 Fig. 2 shows the CV spectrum of the spectra shown in Rg. 1 . 

Fig. 3 shows energy spectra ot>tained by taking the ratios of each measured energy value to the measured energy 
value at wavelength 6000 cm'\ 

Fig. 4 shows ttie CV spectrum of the spectra shown in Rg. 3. 
Fig. 5 shows energy spectra obtained by taking ratios in subdomains. 
so Rg. 6 shows the CV spectrum of the spectra shown in Fig. 5. 

Rg. 7 shows energy specti^a obtained by first-order differentiation of the energy spectra shown in Fig. 1. 
Rg. 8 shows the CV spectrum of the spectra shown in Rg. 7. 

Rg. 9 shows spectra obtained by second-order differentiation of the energy spectra shown in Fig. 1 . 
Fig. 10 shows the CV spectrum of ttie spectra shown in Fig. 9. 
55 Fig. 11 shows energy spectra obtained by ttie low-frequency filtering Fourier transform. 
Fig. 12 shows the CV spectrum of ttie specti'a shown in Fig. 11. 
Fig. 13 shows results of quantitative analysis In the case of no processing. 
Fig. 14 shows results of quantitative analysis in the case of taking ratios for single wavelength. 
Rg. 1 5 shows results of quantitative analysis in the case of taking ratios in subdomains. 
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Rg. 16 shows results of quantitative analysis in the case of first-order differentiation. 

Fig. 17 shows results of quantitative analysis in the case of second-order differentiation. 

Fig. 1 8 shows results of quantitative analysis in the case of performing the low-frequency filtering Fourier transform. 

5 DETAILED DESCRIPTION OF THE PREFERRE D EMBODIMENTS 

The preferred embodiments according to the present invention will be described below with reference to the 
attached drawings. 

10 FIRST EMBODIMENT 

We filled a quartz cell with commeraal milk. Then, using an infrared spectroscope manufactured by Partdn-Elmer 
Co., we measured the energy spectra of transmitted light 1 1 times every 5 minutes on the domain of wavelengths from 
4000 cm-^ to 8000 cm V Then we perfonned, for the measurements of the transmitted energy spectrum, (1) no process- 
75 ing. (2) processing of obtaining the ratio of each measured energy value to the measured energy value at a single wave- 
length. (3) processing of obtaining the ratio of each measured energy value to the measured energy value at a 
predetermined wavelength within each subdomain, (4) first-order differentiation. (5) second-order differentiation, and 
(6) the Fourier transform. Further, we calculated the CV values 100 x standard deviation/mean) of the energy spec- 
trum of transmitted light at each selected wavelength in the wavelength domain of 4000 cm'^ to 8000 cm'\ The spec- 
ie trum obtained by plotting CV values at every selected wavelength Is called a CV spectrum in the following. 

(1) No processing: 

The 1 1 energy spectra of transmitted light measured every 5 minutes of the commercial milk are shown in Fig. 1 . 
25 Their CV spectrum is shown in Fig. 2. 

(2) Ratio calculation for single wavelength: 

The ratio spectra of the 1 1 energy spectra of transmitted light shown in Fig. 1 obtained by calculating the ratio of 
30 each measured energy value to the measured energy value at wavelength 6000 cm"^ are shown in Fig. 3. Their CV 
spectrum is shown in Fig. 4. 

(3) Ratio calculation in subdomains: 

35 We divided the wavelength domain of 4000 cm'^ to 8000 cm'^ of the energy spectrum of transmitted light of the 1 1 
commercial milk measured every 5 minutes into the following 8 wavelength subdomains: 
4000 cm"^ to 4500 cm"^ (4250 cm"^). 
4500 cm"^ to 5000 cm"^ (4750 cm*^). 
5000 cm*^ to 5500 cm"^ (5250 cm'^). 
40 5500 cm-^ to 6000 cm'^ (5750 cm'^). 

6000 cm'^ to 6500 cm"^ (6250 cm*^). 
6500 cm '' to 7000 cm'^ (6750 cm'^). 
7000 cm-^ to 7500 cm*^ (7250 cm"^), 
7500 cm-^ to 8000 cm*^ (7750 cm*^). 
45 Then, we dvided each energy spectrum of transmitted light belonging to each subdomain by the energy value at the 
wavelength shown wittiin the parentheses to obtain the ratio spectra on the whole domain as shown in Fig. 5, Their CV 
spectrum is shown in Fig. 6. 

(4) Rrst-order differentiation: 

50 

The first-order differential spectra obtained from the 1 1 energy spectra of transmitted light nteasured every 5 min- 
utes of ttie commercial milk shown in Fig. 1 are shown in Fig. 8. Their CV spectrum is shown in Fig. 8. 

(Si Second-order differentiation: 

55 

The seoond-order differential spectra obtained from the 1 1 energy spectra of transmitted light measured every 5 
minutes of the commercial milk shown in Fig. 1 are shown in Rg. 9. Their CV spectrum is shown in Rg. 10. 
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(6) The Fourier transform: 

The spectra obtained by the Fourier transform from the 11 energy spectra of transmitted light measured every 5 
minutes of the commercial milk shown in Fig. 1 are shown in Rg. 11. Their CV spectrum Is shown in Rg. 12. 

If we compare the CV spectrum of Fig. 6 with the CV spectra of Fig. 2, 4, 10, and 12. then the CV spectral values 
obtained by ratio calculation in subdomains are smaller, as a whole, than the CV spectral values obtained by no 
processing, ratio calculation for single wavelength, first-order differentiation, second-order differentiation, and the Fou- 
rier transform. Therefore, it is seen that the energy spectra obtained by ratio calculation in subdomains are compara- 
tively stable. 

SECOND EMBODIMENT. 



/5 



20 



25 



30 



35 



40 



45 



SO 



After an examined subject drank a sugar-loaded aqueous test solution (nxxJel name: Trelan 75) manufactured by 
Shimizu Pharmaceutica) Co.. we inadiated a finger of the subject with near-infrared ligtit to measure energy spectra of 
reflected hght eve* y 12 minutes using near-infrared spectroscope (model name: System 2000) manufactured by Parkin- 
Elmet Co. At the same time of the measurement of energy spectra of reflected light, we measured the blood sugar of 
the subiect with a glucose monitor (model name: GM-1320) manufactured by Kyoto Daiichi Kagaku Co.. Ltd.. We con- 
ducted the SMTxiltaneous measurement 44 times in total. Then we performed, for the measurements of the reflected 
energy spectra obtained with the near-infrared spectrosoope. (1) no processing. (2) processing of obtaining the ratk> of 
each nneasured energy value to the measured energy value at a single wavelength, (3) processing of obtaining the ratio 
of each measured energy value to the measured energy value at a predetermined wavelength within each subdomain. 
(4) first-order diflereniation. (5) second-order differentiation, and (6) the Fourier transform. Here, in the atx>ve (2) ratto 
calculation for single Mvelength. we divided each of the following energy values outskje the parentheses by the energy 
value at wavelength 6000 cm'V In the above (3) ratio cak;ulatlon in subdomains. we divided each of the folk>wing meas- 
ured energy values outside the parentheses by the energy value at the wavelength shown within the parentheses. 



4250 cm 
4450 cm* 
4550 cm* 
4650 cm* 
4800 cm 
5500 cm 
5700 cm* 
5900 cm 
6100 cm' 
6400 cm- 
6700 cm* 
6900 cm' 
7200 cm' 
7400 cm' 



(4300 cm*M. 4350 cm'^ (4400 cm'^). 
(4700 cm"^). 4500 cm'^ (4700 cm'^). 
(4700 cm^). 4600 cm*^ (4700 cm'^). 
(4700 cm'M. 4750 cm'^ (4700 cm"^). 
(4700 cm'^), 4850 cm"^ (4700 cm'^). 
(6300 cm*^). 5600 cm'^ (6300 cm"^). 
(6300 cm-^). 5800 cm'^ (6300 cm"''). 
(6300 cm*'). 6000 cm'^ (6300 cm"^). 
(6300 cm'^). 6200 cm'^ (6300 cm'^). 
(6300 cm*^). 6500 cm"^ (6300 cm*^). 
(7100 cm-^). 6800 cm'^ (7100 cm'^). 
(7100 cm*^). 7000 cm*^ (7100 cm'^). 
(7100 cm'^). 7300 cm'^ (7100 cm"*). 
(7100 cm"^). 7500 cm"^ (7100 cm"^). 



Of the data obtained by the above simultaneous measurement conducted 44 times, we used the 22 sets of data 
obtained at even-numbered times for making the analytical curve. We used the 22 sets of data obtained at odd-num- 
bered times for estimating the concernration. In particular, we used the values of ak)sorbance at the 28 wavelengths 
shown outside the parentheses above, and performed quantitative analysis by processing these values of absort)ance 
based on tfie PUS multivariate analysis. 

The results of the quantitative analysis based on (1) no processing. (2) processing of obtaining the ratio of each 
measured energy value to the measured energy value at a single wavelength. (3) processing of obtaining the ratio of 
each measured energy value to the measured energy value at each predetermined wavelength within each sutxiomain, 
(4) first-order differentiation. (5) second-order differentiatk>n. and (6) the Fourier transform are respectively shown in 
Rgs. 13. 14. 15. 16. 17. and 18. The corelation coefficients and estimate errors SEP for the concentration are shown 
in the following Table 1 . 

Here, the estimated error SEP for concentration is defined by the following equation (24). 



55 



SEP. 



2 h,i,ret'^ h,i,obs)^ 



(24) 
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where 

Qi i ref- actual measurement of ooncentration of component / at Mh measurement 
^,hi.obs- estimated ooncentration of component / at h\h measurement. 
n: data size for estimation of concentration. 

5 

Table 1 



Results of Divided Quantitative Analysis 




Least square regression 
equation 


Correlation coefficient R 


Estimate error SEP 
[mgftJQ for concentration 


No processing 


ysr -13.5231 + 1. 08041 X 


0.923082 


14.5359 


Single-wavelength normalize 


y = -10.691 9+1. 05423X 


0.919328 


14.6042 


Subdomain normalize 


y- -19.4918+ 1.11672X 


0.947856 


12.7113 


First-order differentiation 


y = 35.7489 + 0.657103x 


0.669839 


28.3899 


Second-order differentiation 


y = 20.5183 + 0.998778x 


0.882189 


26.3624 


Fourier transform 


y = -12.8506+ 1.07155X 


0.93135 


13.6366 



If we compare Rg. 15 with Figs. 13. 14. 16. 17. and 18. then we see that we can obtain blood sugar estimates 
closer to actual measurements of blood sugar by taking the ratios of energy spectra in subdomains. so that quantitative 
25 analysis of concentration has been improved. Further, from Table 1 , we see that tiie conelation coefficient of between 
blood sugar estin^ted and observed is greater taking tite ratios of energy spectra in subdomains than each con'ela- 
tlon coefficient of t>lood sugar estimates by (1) no processing. (2) ratio calculation for single wavelength, (4) first-order 
differentiation. (5) second-order differentiation, and (6) the Fourier transform. Further, the estimates error tor concen- 
tration SEP is minimum for ratio calculation in subdomains. 

30 

THIRD EMBODIMENT 

We investigated the fluctuations of the surface transmittance p(X.. Q-j due to changes of the incidence angle 0^. when 
light is made incident from air on water, and compared the case In which ratio calculation was performed with the case 
35 in which ratio calculation was not performed. 

We obtained the relative change anxHjnt Ap(X. e|)/p(X.6|) in surface transmittance p(X. Bj due to the change 
amount AO j s e ^ - e of the incidence angle % when ratios were not taken, by the following equation (2^. 

Ap(A.. e,)/p(x.e,) »{p(x, Oj - p(x. e|^),)/p(x. e,,),. (25) 

40 

On ttie other hand, we cbtained the relative change amount Ap,g(X, 6 |)/PN(A..e |) in surface transmittance PfAK Oj) due 
to the change amount Ae , « 6 j ^ - 6 1^ of the incidence angle % when ratios were taken, by the following equation (26). 

apn(a.. eiypN(x.e,) » ip{x. e,.i),/p(Xp e,i), - p(x. e^)^ipix^ e,^),) /[p(x, e, ^),/p(Xr. e, i),). (26) 

45 
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The calculated results are shown in the following Tables 2 and 3. 



Table 2 



10 



IS 



20 



Waveleng/i 
(fim) 






6j(rad) 




1.256 


1.00027 


1.3210 


0 


0.980910 


1.256 


1.00027 


1.3210 


5 


0.981035 


1.256 


1.00027 


1.3210 


10 


0.981416 


1.256 


1.00027 


1.3210 


15 


0.982063 


1.256 


1.00027 


1.3210 


20 


0.982997 


1.256 


1.00027 


1.3210 


25 


0,984242 


0.6563 


1.00028 


1.3311 


0 


0.979864 


0.6563 


1.00028 


1.3311 


5 


0.979996 


0.6563 


1.00028 


1.3311 


10 


0.980396 


0.6563 


1.00028 


1.3311 


15 


0.981076 


0.6563 


1.00028 


1.3311 


20 


0.982056 


0.6563 


1.00028 


1.3311 


25 


0.983363 



25 

Here. 

n^: index of refraction of air. cited from Rika Nenpyo, National Astronomical Observatory (Maruzen Co.), 
02'. index of refraction of water, calculated from the following equation (27) desaibed in Rika Nenpyo. 

30 

{n^ - 1)x10® = 6432.8 + 2949810/(146 - VX^) + 25540/(41 - 1/X^). (27) 



35 

Table 3 



(rad) 


^1.2 

(rad) 


Ap(x. ei)/p(x. 

ej)(rad-i) 


^{^* ei)/PN(>.. 
OiXrad-^) 


0 


5 


1.2743x10-'* 


7.2789x10-^ 


0 


10 


5.1584x10"^ 


2.7070x10'^ 


10 


15 


6.5925x10-^ 


3.4321x10-^ 


20 


25 


1.2665x10'^ 


6.2460x10-5 



As seen from the above Table 3. the relative change amount Apfg(X, 6 |)/PN(X.e ,) in surface transmittance pt^iK 
6i) is smaller than the relative change amount Ap(X. e |)/p(X.e |) in surface transmittance p(X. OJ. Therefore, changes 
so in the surface transmittance p(X. 0) due to changes of the incidence angle 6j are reduced by performing ratio calculation. 
That is. the fluctuations of transmittance due to changes off the incidence angle of light with the measured object can 
be reduced. 

Although the present invention has been fully described in connection with the preferred enrbodiments thereof with 
reference to the accompanying drawings, it is to be noted that various changes and modifications are apparent to those 
55 skilled in the art. Such changes and nxxiifications are to be understood as included within the scope of the present 
invention as defined by the appended claims unless they depart therefrom. 
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Claims 

1. A method of stabilizing spectra in spectrometry that performs quantitative analysis of a specific component in an 
object to be measured by irradiating light to said object and by measuring one of the spectra selected from the 

5 energy spectrum of transmitted light, the energy spectrum of scattered light and the energy spectrum of reflected 
light, comprising steps of 

measuring an energy spectrum of light transmitted through or reflected from said object; 
dividing said energy spectrum into a plurality of wavelength domains, thereby obtaining a plurality of partial 
10 energy spectra; 

normalizing each of said plurality of partial energy spectra within each wavelength domain using an energy 
measured at a predetermined wavelength contained in said each wavelength domain; and 
performing said quantative analysis of a specified component in said object using said plurality of partial 
energy spectra having been normalized. 

IS 

2. The method of stabilizing spectra in spectrometry according to claim 1 wherein said quantative analysis of a spe- 
cific component is performed by multivariate analysis for said plurality of partial energy spectra having been nor- 
malized. 

20 3. The method of stabilizing spectra in spectrometry according to claim 1 or 2 wherein said predetermined wavelength 
is a center value of each wavelength domain. 

4. The method of stabilizing spectra in spectrometry according to claim 1 or 2 wherein said specific component is glu- 
cose. 

25 
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